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Abstract. The 13 — 16/im region observed by the Infrared 
Space Observatory (ISO) of several ABG stars are discussed. 
We present and analyse spectra of seven M giants which show 
carbon dioxide features. To explain the features of the bands, 
we suggest they originate from two different circumstellar lay- 
ers, one being a warm and high density layer close to the star, 
possibly making the 15/im band optically thick, and the other 
being a large, cold and optically thin layer extending far out in 
the wind. This could explain the difference in temperatures of 
the different bands found in the analysis of the spectra and the 
number of molecules needed for the emission. 

It is demonstrated that in spite of the bands probably not 
being formed in vibrational LTE, the temperatures can be esti- 
mated from the widths of the bands. 



Key words: stars: AGB and post-AGB 
- late-type - Infrared: stars 



circumstellar matter 



1. Introduction 

Asymptotic Giant Branch (AGB) stars are considered to be 
important contributors of heavy elements to the interstellar 
medium through their slow but massive stellar winds, and they 
thereby play an important role in the chem ical evolution of 
galaxies, see e.g. Gustafsson & Ryde (1998). The physics of 
this mass-loss and the processes in winds of AGB stars, such as 
the nucleation of dust grains and the acceleration of the wind, 
are, however, not well understood. In trying to gain some in- 
sight into the physical and chemical conditions of the wind re- 
gion and some understanding of the mass-loss mechanism, the 
Infrared Space Observatory (ISO) has been used to observe the 
outer atmospheres and the circumstellar regions of a sample of 
AGB stars by observing spectra of diagnostic molecular bands. 
As one part in this project we have studied carbon dioxide in 
the oxygen-rich, semi-regular asymptotic giant branch stars R 
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Crt, R Dor, R Cas, RS CrA, TX Cam, Villi Oph, and V656 
Cas. 

A dozen triatomic molecules have been detected in AGB 
stars to date, mainly at radio wave-lengths (cf. Olofsson et al., 
1997). From ground-based telescopes gas-phase CO2 has been 
impossible to detect due to the blocking by the terrestrial at- 
mosphere. Only satellite observations provide an opportunity 
to observe astrophysical CO2. 

The first report on CO2 emission from M giants was given 
in Ryde et al. (1997), where it was suggested that the emission 
feature at 14.98/mi in R Crateris was due to the 00°0 E + - 
01*0 II„ Q-branch rotation- vibration transitions of circumstel- 
lar CO2. Based on more data for six M giants, Justtanont et al. 
(1998) could prove unequivocally that this feature and others 
originate from circumstellar gaseous CO2. On the basis of a 
preliminary modelling of the emission lines, Justtanont et al. 



(1998) and Cami et al. (1998) suggest that the emission origi- 



nates from in a region, in which molecules are formed, which 
extends a few stellar radii out from the star. From a thermal 
model they conclude that the gaseous CO2 has a temperature 
of more than 650K. 

Here, we present new ISO observations of the 13 — 15/im 
region of several M giants. CO2 emission is detected to various 
degrees in these stars. 

2. Observations and reductions 

The spectra of the seven M giants and the C star RV Aqr 
were observed with the Short- Wavelength Spectrometer (SWS, 



de Graauw, et al. 1996) on board ISO (Kessler et al., 1996) 



The spectrometer was used in the grating scan mode (SWS06), 
providing a resolution of R ~ 1900, corresponding to A A ~ 
8 • 10 -3 /im or At; ~ 160km/s. Finally one observation of R 
Dor was performed using the high resolution SWS07 mode of 
the Fabry-Perot interferometer with AA ~ 0.5 • 10~ 3 /im. The 
field of view of the SWS06 observation was 14'' times 27". 

The data were partly reduced at the ISOSDC (ISO Spec- 
troscopy Data Centre) at the Max-Planck Institut fur Extrater- 
restrische Physik in Garching, using the 1997 version of the 
Interactive Analysis package (IA), and in Uppsala using ISAP 
1.5 (ISO Spectral Analysis Package). The observed bands are 
not heavily crowded by cosmic ray hits, but are unfortunately 
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Fig. 1. The seven observed M giants and the carbon star RV Aqr as observed by ISO. The data are scaled in strength to make 
a qualitative comparison possible between the stars. Superimposed are the defringed spectra. The observed CO2 Q-bands lie at 
13.9/im and at 15.0/im 

Table 1. Basic data for M giants with CO2 features, according to Loup et al. (1993) 



M giant 


Type 


Class 


Period 

(days) 


M 

(M Q yr 


- 1 ) 


V c 

(kms" 1 ) 


rco 
(cm] 




Distance 

(pc) 


Date of observation 
& revolution no. 


Feature 


RCrt 


M7III 


SRb 


160 


7.5 ■ 10 


-7 


11 


0.68 


■10 17 


250 


960202 (228) 


em 


R Dor 


M8IIIe 


SRb 


338 


0.7- 10 


-7 


6 


0.21 


■10 17 


40 


970623 (585) 


em 






















and 971228 (773) 




RCas 


M7IIIe 


mira 


430 


11 ■ 10" 


7 


12 


0.80 


■10 17 


200 


961203 (382) 


em (?)& abs 


TXCam 


M8.5 


SRb 


557 


25 ■ 10" 


7 


18 


1.2 • 


10 17 


330 


970823 (646) 


abs 


Villi Oph 


M9 


mira 




60 ■ 10" 


7 


17 


1.9 • 


10 17 


570 


970424 (525) 


abs 


V656 Cas 


M9 


mira 




69 ■ 10" 


7 


14 


2.0 • 


10 17 


950 


970928 (683) 


abs 


RS CrA 








28 ■ 10" 


7 


21 


1.1 ■ 


10 17 


280 


961015 (334) 


abs (?) 



severely affected by fringes due to regular interference patterns, 
which depend on the exact pointing of the satellite at the time 
of observation. Therefore, the fringes can not be reduced away 
directly. A possible way to artificially eliminate the fringes is, 
for example, to divide the signal by the simple harmonics that 
dominate the fringe pattern. This was done by using a defring- 
ing algorithm included in the IA package. Such a procedure is, 



however, somewhat questionable since real features, especially 
molecular lines occurring at frequencies close to those of the 
fringe pattern, are also taken away if not special care is taken. 
The flux calibration leads to a relative uncertainty of less than 
10%. 

The M giants observed are given in Table |l]and the observa- 
tions of these stars are shown in Fig. [j]. The observations are all 
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Table 2. Parameters of the CO2 Q-bands. The band oscillator 
strengths were derived from Rothman et al. (1992) and the ro- 
tational constants are directly given in the same reference 
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Wavelength in fim 



Fig. 2. Details of the defringed spectrum of R Dor. The CO2 
Q-bands at 13.9/im and at 15.0/im are asymmetric, the former 
to the red and the latter to the blue 



shown in the same figure and are therefore shifted and scaled 
arbitrarily in vertical direction to enable a comparison between 
the stars. Superimposed on the original, fringed data are the ar- 
tificially defringed data. The carbon star RV Aqr, supposedly 
without any CO2 but with the same instrumental fringes, is 
plotted to give evidence for the hypothesis that the fringes are 
instrumental and not specific for M stars. The fringe pattern 
shows a clear resemblance among the M stars, and is convinc- 
ingly similar between the M giant spectra and the carbon star 
spectrum. The spectra of R Cas and R Crt start from a little less 
than 13.9/im, which means that the 13.9/im band is included 
but is fairly close to the edge of the detected spectrum. This 
part of the spectrum is sampled only once, whereas the central 
regions are covered twelve times. Thus, the edges are more sen- 
sitive to cosmic hits and may therefore be noisier. The spectrum 
of RS CrA is of poor quality due to the low flux measured. 

The CO2 features are found around 13.9 and 15.0/im. Vis- 
ible emission at 13.9/im is seen in the spectra of R Crt and 
R Dor, and possibly in the spectrum of R Cas. The feature at 
15.0/im is seen as an emission in the same stars. In R Cas there 
is, in addition to the emission feature, an absorption part which 
is broader. In TX Cam, Villi Oph, and V656 Cas and possi- 
bly RS CrA the feature is seen mainly as an absorption. In Fig. 
U the 13.9 and 15.0/im emission bands of R Dor show asym- 
metries, the former to the red and the latter to the blue. This is 
seen in all the stars that show CO2 emission. 

3. The carbon dioxide rotation-vibration bands 

The bending-mode, rotation-vibrational band at 667cm~ 1 or 
15.0/im, is a £ + — H u transition of the linear and symmetric 
molecule CO2, see e.g. Herzberg (1966). The selection rules 
for this perpendicular band, At; 



±1 and A J = 0, ±1, result 
in a Q, P, and R-branch. The Q-branch will appear strongest 
since it will correspond to an overlap of A J = transitions, 
whereas the P- and R-branches will be smeared out and will 



wavelength 


fband 




B' - B'Vcm- 1 


13.87/im 


0.167 


• 10" 6 


-1.06 • 10" 3 


14.97/im 


8.65 ■ 


10 -6 


+1.03 • 10~ 3 


16.18/MT1 


0.163 


■ 10~ 6 


-0.77- 10~ 3 



hardly be detectable at the resolution used in the observations. 
If the rotational constant B, determining the energy separa- 
tions of the rotational levels within a vibrational state, would 
be the same for the upper vibrational state as for the lower one, 
the AJ = transitions would all have the same frequency. 
The rotational constant for the upper vibrational level (B') is 
in the 15/im band case, however, modestly larger than that of 
the lower level (B"), leading to a small increase in frequency 
for high J-values. Therefore, the Q-branch is expected to be 
slightly asymmetric to the blue. For the 13.9 and 16. 2/im bands 
(cf. Fig.^J) B' < B" implying a reversed asymmetry. The width 
of a thermal band reveals the underlying temperature since as 
the temperature increases, higher rotational levels will be pop- 
ulated, resulting in a broader band. The difference of the rota- 
tional constant of the upper and lower levels gives an indica- 
tion of how sensitive the broadening is to the temperature. In 
Table I the band oscillator strength, fband and the difference be- 
tween the rotational constants of the upper and the lower level 
of the transition, B' — B", are given for the three bands dis- 
cussed here. 

The 13.9/im Q-band is a transition from a state where the 
molecule vibrates in the v\-mad& and de-excites to the upper 
level of the 15.0/im band. The Q-band at 16. 2/im is a de- 
excitation of a molecule vibrating in the second level of the 
Vi mode down to the upper level of the 15.0/im band. See Fig. 

I 

The resolution of the measurements (AA ~ 0.008/im) 
will allow the asymmetry of the observed features (FWHM~ 
0.015/im) to be measurable, which is also observed. The ro- 
tational fine structure, i.e. the individual Q-lines will, how- 
ever, not be detectable in the SWS06 mode of the spectrom- 
eter. In the high resolution mode (SWS07, AA - 0.0005/im) 
the Q-lines will show up as separated symmetric lines, since 
the typical separation for the lowest Q-lines is approximately 
0.001/im. 

4. Results and analysis 

Observational evidence points toward a relative simplicity of 
the kinematics, the density and the temperature structures of 
stellar envelopes, with an over-all symmetry prevailing (see e.g. 
Olofsson, 1996), although several AGB stars have been shown 
to have asymmetric winds (Plez & Lambert, 1994 ). Also, note 
the clear evide nce o f clumpy structure at some scale in e.g. 
Olofsson et al. ( 1998 ). The chemical structure is certainly quite 
complex, with a chemistry far out of equilibrium. A standard 
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Table 3. Maximum optical depths, r n 
for a spherical and stationary wind 
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Fig. 3. The vibrational energy-level diagram of some states in 
CO2. A number of infrared transitions are indicated. The labels 
identifying the different levels indicate the number of quanta in 
the three fundamental vibrational modes. The superscript de- 
notes the quantum number associated with vibrational angular 
momentum, I 

circumstellar envelope, with a non-rotating, smooth, and spher- 
ically expanding wind with a constant velocity is nevertheless 
a reasonable starting point for our analysis. 

4.1. Are the bands optically thin? 

In the following we will discuss the optical thickness of the 
environment from which the bands are emitted. The order-of- 
magnitude estimate made here will lead to an upper limit partly 
because the wind velocity will eventually shift the molecules 
out of the thermal broadening of the Q-branch lines due to 
Doppler shifts and thereby lower the estimated optical depth. 
For an estimate of this upper limit of the optical depth, t, of 
the bands, we assume a spherically symmetric wind resulting 
from a constant rate of mass-loss. Within the region where CO2 
is assumed to exist, its abundance relative to hydrogen, e co , 
is set constant. The optical depth is found by integrating the 
cross section per particle (ax, [cm 2 ]) for the absorption times 
the number density of the absorber (N, [cm -3 ]) along the line 
of sight with an impact parameter p. For atomic or molecular 
scattering <j\ = 7re 2 A 2 fb a nd/(nieC 2 <5A), with fband being the 
band oscillator strength. From the continuity equation and the 
assumption of a stationary flow the number density of CO2 in 
all energy states is found. Thus, in cgs units, 



T{p) = 21 • 10 9 e c 



MA 2 f, 



band 2 . \/ Rq ' 

— — arctanf— 

V e 0A p p 



), (1) 



where R is the outer boundary of the emitting sphere and 
fi5/im = 9 ■ 10 -6 . The largest optical depth will be found in 
the 15.0^m band since the lower lying energy state of this 
transition is the vibrational ground state, where all but a few 
molecules will most probably be found. An upper limit of the 
optical depth will therefore be estimated for an impact param- 



M giant 


maximum optical depth, r max 


RCrt 


0.07 


R Dor 


0.01 


RCas 


0.09 


TXCam 


0.13 


Villi Oph 


0.33 


V656 Cas 


0.05 


RS CrA 


0.13 



eter p — 1 R* and assuming that all CO2 is in the ground 
state. An upper limit of the fraction of CO2 (e C02 ) is found 
by adopting a reasonable number of the CO2/CO ratio and as- 
suming that almost all carbon is locked in CO. This is not a 
bad assumption in an oxygen-rich environment since the CO 
molecule has a large binding energy and is the most stable 
molecule in the circumstellar environment. The CO2/CO ra- 
tio is assumed be a few times 10 -4 , in accordance with the 
calculated abundances in the circumstellar chemical models of 



Willacy & Millar (|1997j). The width of the band (SX) is approx- 
imately 0.02/im. 

The results for the stars considered are tabulated in Table 
[I which gives the optical depths for a standard wind extend- 
ing out to around 1000 stellar radii. Obviously, all bands will 
be optically thin in this case. A wind extending out only to 3 
stellar radii, would give around 20% lower optical depths un- 
der the same assumptions. If we consider a warm and quasi- 
static, molecular-forming region with small velocities, such as 
the one suggested by Tsuji et al. (1997), the densities close to 
the star are enhanced, which will lead to higher optical depths. 
A density enhancement close to the star by one or two orders of 
magnitude could make the 15/im band optically thick. An in- 
crease of the density could be achieved if extra ejected matter 
eventually falls back on the star without raising the net matter 
transport outwards. Turbulent velocities will broaden the lines 
and the quasi-static layer will not shift the molecules out of the 
lines. Thus, there is a possibility that the absorption bands ob- 
served are formed in such high density regions. In view of their 
lower oscillator strengths and their excitation, the other bands 
are obviously optically thin also in this picture. 

4.2. Are the bands formed in LTE? 

It is widely believed that most AGB stars pulsate with shocks 
travelling through their circumstellar regions. These dynamic 
regions are often assumed to show departures from LTE condi- 
tions. We will show that it may be questionable to assume that 
the bands are in LTE. 

The highest probability for LTE to prevail, in an expand- 
ing wind model, is close to the star since the density de- 
creases, at stationary outflow, with the square of the distance 
from the centre. From the equation of continuity for a spher- 
ically symmetric wind we get the number density of hydro- 
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gen molecules as a function of the radial distance from the 
star r, as A H2 = Af/(47rw c r 2 m H2 ). For R Crt with M = 
8 • 10" 7 M Q yr- 1 and v c 



M = 7 ■ lO- 8 M yr- 1 and v, 



1 1 km s 1 and for R Dor with 
6kms -1 (Loup et al., 
1993), we get for the number density of H2 for this station- 
ary model, adopting R» = 3 • 10 13 cm: 1 • 10 9 • (i?»/r) 2 cm~ 3 
and 0.2 • 10 9 • (i?*/r) 2 cm~ 3 , respectively. 

Observe, however, that the atmospheres of these stars ex- 
tend far out due to pulsations and shocks, which will invalidate 
the assumption that the wind is in a steady state at small radii. 
Thus, from dynamic models of variable stars by, e.g., Bessel 
et al. (119961) and Hofner & Dorfi (|1997|), it is indeed seen that 



the densities at a few stellar radii from the star can be, in gen- 
eral, an order of magnitude higher than the ones given by the 
equation of continuity. Note also that observations indicate that 
the winds may be inhomogeneous and clumpy (e.g. Olofsson 
et al., 1998), which could increase the densities in the emis- 
sion regions. Also the quasi-static layer discussed above will 
have higher densities, which could be enough for collisions to 
be dominant a few stellar radii out. However, from the sim- 
ple considerations above we will see that we may nevertheless 
conclude that the assumption of statistical equilibrium may be 
needed instead of that of LTE, even a few stellar radii away 
from the star. This must be considered when discussing the ki- 
netic gas temperature and the site of formation of CO2. 

The statistical equilibrium equations, including radiative 
and collisional processes for a set of transitions, give the popu- 
lations of the levels. These will also be given by the Boltzmann 
distribution if the density is sufficiently high in order for col- 
lisions to be dominant over radiation processes. The spectral 
lines are then formed in LTE and the gas is thermalised at the 
kinetic gas temperature. If the collisions are not dominant, it is 
not clear whether LTE is valid. 

An order-of-magnitude estimate of the probability for col- 
lisional excitation is given by C\ u ~ A H2 a v. The cross sec- 
tion, a, can be approximated with the geometrical area of a 
molecule, which is about er gC om = 10~ 16 cm 2 . The velocities 
of the molecules for the temperatures considered are of the 
order of v ~ lkms -1 . For the Q-band at 15.0/^m, with an 
Einstein coefficient for spontaneous decay of A u ] ~ Is -1 , the 
probability of radiative excitation (Bi u • J„, J u being the mean 
intensity of the star) is 0.3s -1 for a typical M giant and ex- 
trapolating back to 1 R*. In order for the collisional processes 
to dominate we thus need A H2 ^> 3 • 10 10 cm -3 . Therefore, it 
seems that the densities in the circumstellar environments, are 
not high enough for the collisions to be completely dominant. 
Radiative processes may therefore be more important than col- 
lisional ones, which may lead to severe departures from LTE. 

However, among the rotational energy levels in the vibra- 
tional ground state the levels are thermally populated. This con- 
clusion can be drawn since the Einstein coefficient for spon- 
taneous decay between the rotational states are several orders 
of magnitude smaller than the coefficients for vibrational de- 
cays. The relative population of the levels within the vibrational 
ground state are therefore given by the Boltzmann distribution 



and the local gas temperature. The population of the rotational 
levels of higher lying vibrational states are affected by radiative 
transitions between the different vibrational states. 

In a situation dominated by radiative processes, the 
strongest transitions involved in a case of statistical equilibrium 
are shown in Fig. |3|. These levels and the transitions between 
them were included in a Monte Carlo simulation of a radiative 
situation to discern to what extent the transitions will eventu- 
ally spread out the populations of the rotational levels within 
higher lying vibrational states. In this situation the rotational 
levels of the ground state are given by the Boltzmann distribu- 
tion and the probabilities of the transitions, i.e. the oscillator 
strengths, of an individual, rotation-vibration transition within 
a band are given by 



rabs 
line 



s. 



J" J' 



2, J" + 1 



A,-. 



2 - S ,_ 



2 - <Va'+A" ' 



(2) 



where A denotes the projection of the electronic orbital 
angular momentum along the axis of the linear molecule and 
Sj"j' ar e the H onl-London factors which can be found in, e.g. 
Bernath ( 1995 ). These can be calculated for all the transitions 
considered. 

The result of this investigation shows that the spreading out 
of the level populations is on a few percent level and hardly 
noticeable. This means that also in an environment dominated 
by radiative processes the relative rotational populations within 
the higher lying vibrational levels will map the Boltzmann 
distribution of the ground vibrational level. This leads to the 
conclusion that the widths of all vibration bands will reflect 
the underlying gas temperature. However, in general, the rela- 
tive strengths or amplitudes of the bands will not be given by 
the thermal distribution for the vibrational levels, if the radi- 
ation processes dominate. The 13.9 and 16.1/im bands origi- 
nate from upper states of about the same energy and they are 
affected similarly by radiation excitation, which may, in any 
case, lead to similar band strengths. 

4.3. The temperature of the bands 

In this section we account for the excitation temperatures of 
the different bands for the stars showing emission. This is done 
by modelling the band widths of a thermally populated CO2 
molecule and fitting its spectrum to the observed bands. The 
HITRAN molecular database ( feothman et al., 1987] ) and the 
high temperature analogue HITEMP database (Rothman et al., 
in prep.) are used to retrieve relevant spectral line data such as 
line positions, level populations and transition probabilities. 

The widths of the bands depend on the population of the ro- 
tational levels within the upper vibrational state and thereby the 
excitation temperature. In a thermal model these are populated 
according to the Boltzmann distribution. 

The observations of the 15.0 and the 13.9/im Q-bands of 
carbon dioxide in R Dor and R Crt are shown in Figs. Q and |5| 
These spectra show the residual emission after the fringes, due 
to the spectrometer, have been removed. Also plotted in the 
figures are the bands of a thermal model for 650K smoothed to 
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Fig. 4. The full line shows the two observed Q-bands of R Cra- 
teris and the dashed line represents a thermal model of 650K 
scaled to the 13.9/im transition 



o 
o 
o 



R Doradus 



13.0 13.5 14.0 14.5 
Wavelength in fj,m 



5.0 



1 5.5 



Fig. 5. As in Fig. Q but for R Doradus. The forests at 14.5 and 
15.5/im are the R and P branches of the 15.0/im transition in 
the model of 650K. These lines are not resolved in the resolu- 
tion of the observation 



the resolution of the observations. The emission predicted by 
the model is plotted with the observed continuum as a base line 
and is scaled to fit the 13.9/im band. 

For R Dor the excitation temperature as estimated from the 
13.9/im band width is T cx = (650 ± 50)K, while the 15.0/im 
band is fitted best with T ox = (150 ± 50)K. The band ratio 
between the 13.9/im and the 15.0/im bands is 1.7. 

Also for R Crt the 13.9/im band suggests T ox = (650 ± 
50)K, whereas the 15.0/im band is better fitted with a thermal 
model with a temperature of T ex ~ (160 ± 50)K, cf. Fig. ^. 
The observed band ratio is for this star 1.0. 

It is also clear from Figs. Q and |5] that the band intensity 
ratios can not be explained by a simple, optically thin, one- 
component LTE-model. Our observations of the 15.0/im band 
of R Dor in the high resolution mode also show the individual 



Fig. 6. The observed 14.98/im-band for R Crt is plotted in 
full line. Arbitrarily scaled thermal models for temperatures of 
125K, 160K, 250K, and 650K are also plotted in the figure. The 
maximum of the model spectra is shifted towards the blue with 
increasing temperature. The relative errors in the observations 
are on the order of a few percent 
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4.5 - 



g 4.0- 
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Z 3.5 - 
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2.5 
1 4.94 




14.95 



14.99 15.00 



14.96 14.97 14.98 
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Fig. 7. The high-resolution Fabry-Perot observations of R Dor. 
These observations are very difficult to reduce due to cosmic 
ray hits, and spurious features are a problem. Also plotted in the 
figure is the spectrum smoothed to the resolution corresponding 
to the other observations 



Q-lines. These observations are unfortunately very difficult to 
reduce due to cosmic hits. However, from the observations, it 
can be concluded that the intensity distribution of the different 
rotational lines within this vibrational band suggest a tempera- 
ture several hundred degrees below 650K, cf. Fig. 0. At 650K 
the maximum line strength lies around the transition Q24, and 
at 150K at Q10. Also, the Fabry-Perot observations smoothed 
to a resolution corresponding to the one used in the other spec- 
tra, can not be fitted with a 650K model, cf. Figs. ^ and ^[ 



N. Ryde et al.: The nature of the circumstellar CO2 emission from M giants 



7 



Table 4. Number of excited CO2 molecules required for the 
measured flux and the total number of molecules for two dif- 
ferent shell sizes. See text 



M giant 


Flux/ 

ltr 11 

erg s _1 cm -2 


Nco 2 / 
10 44 


10 44 
(3R.) 


10 44 

(large shell) 


RCrt 










13.9/im 


3.9 ±0.4 


350 ± 40 


0.6 




15.0/im 


2.9 ±0.3 


12 ±1 


0.6 


300 


R Dor 










13.9/im 


21 ±2 


48 ±5 


0.09 




15.0/im 


9.5 ±1 


0.97 ±0.1 


0.09 


50 


RCas 










15.0/im 


2.3 ±0.3 


5.9 ±0.6 




400 



4.4. The number of emitting molecules 

The number of molecules needed in a certain vibrational state 
in order for an emitting optically thin envelope to radiate a 
given measured flux, F m , can be calculated. With the Ein- 
stein coefficient for spontaneous decay, A u i, the total number 
of molecules in the upper level, 7Y tot u , is given by 



4-7T • distance 2 • F„ 
hi/A T1 i 



(3) 



This number should be compared to the total number of 
available CO2 molecules in all vibrational states, which can be 
estimated from 



N(C0 2 



M 



-(Ro — Ri) e co 



[CQ 2 ] 
[CO] ' 



(4) 



where e co is the abundance of CO relative to hydrogen, M the 
mass-loss rate, v c the terminal wind velocity and R and Ri 
are the outer and inner radii of the emitting shell, respectively. 
As before, we adopt a number ratio of CO2 molecules to CO 
molecules of a few 10~ 4 . 

Using these expressions we can compare the number of ex- 
cited CO2 molecules needed for the measured flux to the total 
amount of CO2. The measured fluxes in the 13.9 and 15.0/itn 
bands are presented in Table [|. The table also gives the num- 
ber of excited molecules required for the measured flux as well 
as the total number of molecules in an emitting shell, 3 stel- 
lar radii wide, and a large emitting envelope. The extension of 
the CO2 shells for R Dor, R Crt and R Cas used for the results 
in the last column in Table f|, are found by scaling the results 
for TX Cam from the models of Willacy & Millar. This column 
therefore gives the number of carbon dioxide molecules in a 
large emitting shell extending from 50 to 1000 R* for R Crt 
and R Cas and from 20 to 300 R* for R Dor. 

5. Discussion 

The excitation temperatures derived from the different bands 
are not in agreement with each other. We therefore suggest a 
two component model which could explain the spectra. 



In the first component, corresponding to the region sug- 
gested by Justtanont et al. (1998), the CO2 exists in a dense 



quasi-static layer close to the star, the layer extending out 
to only a few stellar radii, resembling the warm, molecular- 



forming region suggested by Tsuji et al. (1997). CO2 is sup- 
posed to have a temperature of the order of 650K and high 
densities. This would give rise to the 13.9/itn band and the 
absorption at 15/im, a band which is now optically thick and 
which will appear differently depending on the physical condi- 
tions. As seen in Table the number of molecules needed for 
the observed emission is too high compared to the number cal- 
culated from a spherically symmetric wind at a constant mass- 
loss rate. It is impossible to produce these fluxes from a small 
region close to the star extending out to, e.g., 3 stellar radii 
with the assumptions above. The number of excited molecules 
needed exceeds the expected number for a thermally populated 
level by several orders of magnitude. 

More reasonable numbers will result from a high density 
layer with a density of more than Ah 2 = 10 10 cm~ 3 , only ex- 
tending out to a few stellar radii. A density enhancement of 
several orders-of-magnitude above the classical wind density, 
is required. In this case, the 15/im band will be formed in con- 
ditions in which collision are c lose to being dominant over ra- 



diative processes, cf. Sect. 4.2 



The second component, which has been argued for by Ryde 



et al. (1997), would be a region extending far out from the star. 
The colder 15/im emission is then supposed to originate from 
this tenuous region. The corresponding cold 13.9/im emission 
would be hardly detectable in the 650K emission. Thus, as is 
seen in Table assuming a much larger shell of CO2, extend- 
ing out to several hundred stellar radii (finally limited by the 
photo dissociation by the interstellar UV field) the required de- 
gree of excitation in the upper level of the 15.0/im band would 
be some percent for all our stars, which is not unreasonable. 
The number of stellar IR photons per square centimetre and 
second is typically two orders of magnitude greater than the 
number of photons in the measured lines, which is reassuring. 

The temperature and number of molecules required from 
the observed 15/im emission features are thus consistent with 
the location of the carbon dioxide in a circumstellar shell, ex- 
tending typically out to 10 3 stellar radii. This is also consis- 
tent with the chemical models of Willacy & Millar (1997) for 
oxygen-rich circumstellar envelopes, in which CO2 is formed 
in a region from around 100 to 1000 stellar radii. Abundances 
of other molecules in their models show good agreement with 
observations, which is reassuring. In these models carbon is 
provided by CO and CH4, and most carbon-bearing molecules 
are due to the breakdown of CH4. An important reaction for 
the present discussion is the one where some CO is destroyed 
at small radii by the reaction CO+OH— >C02+H. Further out, 
all molecules will eventually be destroyed by photo-dialysis by 
the interstellar ultraviolet field. 

R Cas is a special case (Fig. [j]). Here, both absorption and 
emission are clearly seen. With the emission from the cold, 
second component filling in the absorption at 15/im in the red 
wing, the width of the absorption is consistent with 650K. The 
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residual emission in the 15/im band, after the absorption has 
been subtracted, is (2.3 ± 0.3) • 10 _11 erg s _1 cm -2 , which re- 
quires (5.9 ±0.6) • 10 44 molecules. This is a reasonable number 
of excited molecules compared to the total number expected in 
a large shell, cf. Table [|. 

As seen from Table [l], where the parameters of the stars 
discussed are given, the only parameter that is clearly different 
between the M giants with and without CO2 emission is the 
mass-loss rate. In the latter stars the high mass-loss rate seems 
to prevent the formation of the CO2 emission. This finding was 



solved Q-band. However, in analysing our seven ISO spectra of 
M giants, we find that the gas temperatures may be estimated 
from the widths of the bands even if the bands are not in vi- 



brational LTE. We note that Onaka et al. (1998) show that all 



also hinted at by Justtanont et al. (1998). They found that the 
equivalent widths of the CO2 feature are well and positively 
correlated with the 13/im dust feature. They indicate that the 
higher the mass-loss rate the smaller is the 13/im feature and 
thereby the less pronounced are CO2 features. This is consis- 
tent with our findings. Whether the correlation between mass- 
loss rate and the CO2 features has a simple physical origin is 
still uncertain. It seems, however, tempting to suggest that the 
stars with high mass-loss rates more easily accelerate winds 
from the stellar surface to velocities above the velocity of es- 
cape while stars with lower mass-loss rates produce less ac- 
celeration. For these stars much of the material levitated above 
the photosphere is stopped by gravity or shocks close to the star 
where it forms a relatively dense warm layer, where CO2 and 
dust may form and radiate. 

The relative strengths of the 15.0/im band emission and the 
13.9/im band emission are different in R Crt and R Dor. This 
difference in rel ative strengths is also seen in the six stars of 
Justtanont el al. ( 1998 ). This would be a natural consequence of 
the two component model. In an optically thin, one component 
model all molecules de-exciting from the 10°0 and 02°0 energy 
states will finally have to de-excite via the 15.0/im transition. 
This does not seem to occur according to the observations. A 
two component model with one of the bands being optically 
thick may solve this. It could also be explained if the resonance 
band is on the verge of becoming optically thick, thereby being 
weakened by self absorption. But in this case one would sooner 
expect a broader band which is not seen in the observations. 
Our high resolution observations also do not verify this. 

6. Conclusions 

The temperatures deduced appear not to be consistent with a 
model of one single temperature component. Therefore, we 
have suggested that the bands originate from two different com- 
ponents: One would be a warm and high density layer close to 
the star, possibly making the 15/im band optically thick, an- 
other would consist of a large, cold and optically thin CO2 
region extending out in the wind to hundreds of stellar radii. 
These two components could explain the differences in tem- 
peratures and in the number of molecules found from the ob- 
servations. 

The circumstellar CO2 emission bands observed in M gi- 
ants are probably formed out of LTE and caution should be 
exercised when discussing the gas temperature and the site of 
formation of CO2 from observations of the width of an unre- 



bands correlate with variability phase. This could be explained 
by differences in the radiation field if the population is deter- 
mined by radiative processes. 

Since the CO2 molecules may be dominated by radiative 
pumping by stellar flux via one or some of the many vibrational 
levels, a detailed treatment of CO2 in the radiation field and 
wind of the star seems important for modelling the emission 
and verifying the hypothesis of its origin. 
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